First order metamagnetic transition in Ho 2 Ti 2 07 observed by vibrating coil 
magnetometry at milli-Kelvin temperatures 
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We report vibrating coil magnetometry ol the spin ice system Ho2Ti207 down to ~ 0.04 K for 
magnetic fields up to 5 T applied parallel to the [111] axis. History dependent behavior emerges 
below Tq ~ 0.6 K near zero magnetic field, in common with other spin ice compounds. In large 
magnetic fields we observe a magnetization plateau followed by a hysteretic metamagnetic transition. 
The temperature dependence of the coercive fields as well as the susceptibility calculated from the 
magnetization identify the metamagnetic transition as a line of first order transitions terminating 
in a critical endpoint at ~ 0.37 K, B m ~ 1.5 T. The metamagnetic transition in Ho2Ti20y is 
strongly reminiscent of that observed in Dy2Ti207, suggestive of a general feature of the spin ices. 

PACS numbers: 75.30.Kz; 75.60.Ej; 75.40.Cx; 75.40.Gb 



Metamagnetism (MMT), where a paramagnet under- 
goes a first order phase transition to a ferromagnetic 
(FM) state in high magnetic fields with a discontinuous 
jump in the magnetisation, is a pervasive phenomenon in 
systems based on rare earth or transition metals. How- 
ever, despite striking similarities in the magnetization 
of a wide range of MMTs their microscopic origin may 
be radically different (see e.g. PUS]). A prominent ex- 
ample of such a MMT has been reported recently for 
Dy2Ti 2 07, where the transition is preceded by a plateau 
in the magnetization and a critical end-point is located at 
T m ~ 0.36 K and B m ~ 0.9 T for fields strictly parallel to 
a (111) axis [3 [8]. It has been argued that the MMT in 
Dy2Ti2C>7 reflects directly the nature of the spin excita- 
tions from the zero- field spin state [9] . The magnetic ions 
in Dy2Ti2C>7 reside on the vertices of a pyrochlore lat- 
tice of corner sharing tetrahedra in the presence of strong 
local (111) crystalline anisotropy and effective ferromag- 
netic interactions, leading to geometric frustration. The 
ground state is characterized by a residual entropy quan- 
titatively comparable to the value of water ice [T0lU2| . 
This reflects spin disorder at low temperature such that 
two spins are constrained to point outward and two spins 
in towards the center of a tetrahedron (TTJ [13] . 

The MMT in Dy 2 Ti 2 07 has been explained in terms of 
an entropy reduction which takes place in two steps |14j . 
First, the system partially magnetises retaining the two- 
in, two-out state, in which one of the four spins on each 
tetrahedron has a component of the moment antiparal- 
lel to the field. Second, as the field increases further, 
the nearest-neighbor spin ice model predicts an ice-rule 
breaking spin flip to the three-in, one-out (one-in, three 
out) state [13 HE] ■ Recent theoretical work suggests that 
the spin-flips related to the second step may be viewed 
as emergent magnetic monopoles that condense at the 
MMT [HI E] • This scenario was subsequently found to 



be consistent with the entropy reduction inferred from 
the magnetocaloric effect [U HH] , as well as the evolution 
of the spin relaxation time in the ac susceptibility [T9|l20j. 
at least above 1 K. The strong temperature dependence 
of the specific heat below 1 K [3T] , the heat transport 
[22] and finally magnetization avalanches in low mag- 
netic fields [33] are also thought to be consistent with 
magnetic monopoles. On a microscopic level, evidence 
for magnetic monopoles in Dy2Ti 2 C>7 has been inferred 
from the observation of lines of reversed spins between 
monopole pairs (cf. Dirac strings) using neutron scatter- 
ing [3T]. 

A second candidate for emergent magnetic monopoles 
is rIo 2 Ti 2 07. The experimental situation in Ho2Ti 2 07 
is, however, much less clear. Neutron scattering at low 
temperatures reveals pinch-points in the structure factor 
consistent with power law correlation functions [3H 135] . 
However, neutron spin echo and neutron backscatter- 
ing experiments [3B-- 28 suggest intrinsic spin relaxation 
times much faster than for Dy 2 Ti 2 07, raising the ques- 
tion how this may be reconciled in terms of magnetic 
monopoles. Further, the specific heat displays large nu- 
clear hyperfme contributions [29], which manifest them- 
selves as a Schottky anomaly, complicating a comparison 
with the predictions for spin ice behaviour. The magne- 
tization of Ho2Ti2C>7 for a field along [111] reported so 
far down to 0.5 K [301 [3T] showed a non-linear increase 
around 2T reminiscent of Dy 2 Ti 2 07 (distinct from dis- 
cussions of a liquid-gas transition for fields along [100] 
[32]). In general, the magnetic phase diagram of Ho- 
based compounds may show strong effects below 0.5 K 
due to hyperfme interactions as e.g., for the transverse 
field Ising magnet LiHoF4 [33] . Hence, detailed measure- 
ments well below 0.5 K are needed to explore the putative 
equivalence between Ho 2 Ti207 and Dy 2 Ti 2 07 quantita- 
tively. 
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The observation of history dependencies and dynamics 
on long timescales in Dy2Ti207 and Ho2Ti207 in pre- 
vious studies imposes strong experimental constraints. 
First, tiny sample movements in the applied magnetic 
field must be avoided; they might change the magnetic 
state. Second, the sample must be rigidly anchored ther- 
mally, since changes of the magnetization may cause large 
associated entropy release and uncontrolled local heating 
effects. Third, as the magnetic properties are sensitive to 
the precise field value and orientation [341 135] a uniform 
applied field is essential. 

All of these requirements are met by the vibrating coil 
magnetometer (VCM) [36l [37] we used to measure the 
magnetization of Ho2Ti 2 07. Data reported in the fol- 
lowing between 0.04 K and ~ f.8K correspond to the 
properties for field parallel to [If 1] within a few tenths 
of a degree. For larger misalignments (~ ±2°) we found 
that the features of interest tend to broaden and shift, 
with additional hysteretic features suggesting complexi- 
ties beyond the scope of our study. All temperature de- 
pendent data were recorded while continuously heating 
at a rate of 5 mK/min, where zero-field cooled properties 
and field-cooled properties are distinguished. Likewise, 
all measurements as a function of magnetic field followed 
one of two well defined protocols, denoted as (A) and (B). 
Details of the temperature dependent measurements as 
well as the two protocols are presented in the supplemen- 
tary material. 

The Ho2Ti 2 07 single crystal studied was grown by 
means of optical float zoning at LANL in an argon tmo- 
sphere at 5mm/hr. Single crystal neutron diffraction at 
HEIDI (FRM II) establish a room temperature lattice 
constant of a=10.13(2) A consistent with previous stud- 
ies [35]. A psi scan of the (555) reflection confirms the 
homogeneity and high quality of the sample. Demag- 
netising fields were accounted for by approximating the 
disc shaped single crystal as an ellipsoid (7.3 x 4.8 x 1.2 
mm 3 ) with a demagnetising factor N=0.75 ([111] direc- 
tion perpendicular to the plane) [35J 00] ■ The edges of 
the sample were wedge shaped, which may result in a dis- 
tribution of internal fields over a small volume fraction. 
This does not affect the main conclusions reported here. 

Shown in Fig.[T](a) is the temperature dependence of 
the magnetization of Ho2Ti207 for B = 0.01 T along 
[111]. With decreasing temperature the low-field magne- 
tization increases gradually, characteristic of a paramag- 
netic state. Below a temperature Tq ~ 0.6 K the zfc/fh 
and fc/fh data begin to show pronounced differences. 
While the former rapidly decreases to a very low value, 
the latter remain constant. This behaviour is strongly 
reminiscent of other spin ices, where Tq is roughly 0.6- 

0. 75 K for A 2 B 2 7 (A=Dy,Ho; B=Ti, Sn) [30l SB [42], 

1. e., Tq is not proportional to the rare earth exchange 
coupling J |43j . The microscopic origin of the history 
dependence shares many features of magnetic blocking, 
but is not understood. For instance, Monte Carlo simu- 
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FIG. 1: (Color online) Characteristic features of the mag- 
netization of spin freezing, hysteresis and a metamagnetic 
transition ending in a critical point in Ho2Ti207. (a) 
Temperature dependence of the magnetization of Ho2Ti207 
for [111] in a small applied field of 0.01 T. (b) Temperature 
dependence of the coercive fields with respect to zero field, 
(c) Temperature dependence of the coercive fields at the 
metamagnetic transition, (d) Temperature dependence of 
the peak value of the inverse susceptibility with decreasing 
temperature, approaching the metamagnetic transition. 



lations predict a fst order transition at 0.18K [H] which 
has never been observed experimentally. Interestingly, 
however, an exponential slowing down of the spin relax- 
ation has been reported in solid paramagnets which pro- 
ceeds via energy levels caused by crystalline or hyperfine 
splitting of the ground state of the ion [4"5] . 

As a function of magnetic field the history dependence 
below T * is connected with strong hysteresis with respect 
to B = 0, followed by a metamagnetic increase at a field 
B^ ~ 1.5 T, which becomes distinctly hysteretic at low 
temperatures. This is illustrated in Fig. [2] where typical 
data are shown as a function of internal magnetic field, 
-Bint, given by [IqH+M. Data at f .8 K shown in Fig.[2](a) 
are in excellent agreement with Refs. [301 131] ■ The vari- 
ation of the magnetization near B = and around f.5T 
becomes much steeper below T * = 0.7 K as shown in 
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FIG. 2: (Color online) Magnetic field dependence of the low 
temperature magnetization of Ho2Ti207. With decreasing 
temperature, hysteresis emerges at B = and around 1.5 T. 
The susceptibility calculated from the magnetization is shown 
on the right hand side of each panel. At the metamagnetic 
transition it increases strongly with decreasing temperature. 
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FIG. 3: (Color online) magnetization of Ho2Ti207 as a 
function of internal magnetic field after correction of the 
demagnetising fields, following field sweeps of protocol (B). 
For T < Tq the initial change of the magnetization is 
zero within experimental sensitivity and followed by a very 
pronounced increase with increasing field. (a) Magnetic 
field dependence in sweeps of type Bl. (b) through (d) 
magnetization in field cycles of sweep types Bl, B2 and 
B3 (e) magnetization in field cycles of sweep types A4 and A5. 



Fig.[2](b). There are several regimes: essentially no hys- 
teresis is observed at 0.7 K within the accuracy of our 
set-up. At T = 0.4 K < T*, shown in Fig.[2](c), size- 
able hysteresis may be seen with respect to zero field, 
while the magnetization at the second step rises non- 
hysteretically, but more steeply than at 0.7 K. Finally, as 
shown in Fig.[2](d), hysteresis exists with respect to both 
B = and B m = 1.5 T for 0.1 K. 

In order to track the width of the hysteresis loop wc de- 
fine coercive fields Bq , Bq and (see Fig.[2](d)). 
Shown in Fig.[l](b) are the coercive ields Bq , Bq , which 
increase strongly below Tq with decreasing temperature 
and appear to level off around ±0.2 T as T — > 0. In con- 
trast, the hysteretic behaviour at high fields appears at 
T* n = 0.37 K, well below T * (Fig.Qc)). 

No conventional long range thermodynamic ordering 
transition is reported as a function of temperature at 
B = 0, despite the appearance of strong hysteresis. 
Specifically, there is no development of magnetic Bragg 
peaks. Rather, the notion of strong magnetic blocking 
is beautifully illustrated by the initial field dependence 



of the magnetization for T < T * after zfc. Shown in 
Fig.|3](a) is the magnetization in the first field sweep af- 
ter zero-field cooling (referred to as Bl). Up to field 
values exceeding the coercive fields seen in field sweeps 
(B2), (B3) the magnetization is unchanged and trapped 
in the zfc state. The peculiar field dependence, specif- 
ically the negative slope of the magnetization at in- 
termediate fields is reminiscent of Dy2Ti207 (Fig.[3](b) 
through (e)), where it has been interpreted as evidence 
of monopole avalanches |23j . However, similar features 
have also been reported in mesoscopic spin systems [46] . 
In view of the enormous sensitivity to the precise orienta- 
tion of the sample and the large demagnetising factor, as 
well as a small distribution of internal fields at its fringes, 
a detailed account of possible magnetization avalanches 
is beyond the scope of our study. 

Two arguments establish that the hysteresis at B m 
is connected with a thermodynamic phase transition, in 
contrast to the behaviour in zero field. First, the in- 
verse susceptibility at B m calculated from the magneti- 
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zation essentially displays a Curie dependence and van- 
ishes within experimental accuracy at = 0.37(0.01) K 
and B* m = 1.52(0.01) T as shown in Fig.[l](d). This is 
characteristic of a critical point at (T^-B^J. Moreover, 
the hysteresis observed for T < T* n provides evidence 
that this critical point is located at the end of a line of 
first order transitions. Interestingly, T* n is found to be 
very close to the value for Dy 2 Ti 2 07 [7 . By contrast, 
theory predicts scaling of B* n with the effective exchange 
coupling Jeff PUSH- Given J off = 5D/3 + J/3, where D 
is the dipolar coupling between the rare earth moments, 
and using J e ff/fcB = 1.1 K for Dy2Ti207 and 1.83 K for 
Ho2Ti207 pj], the phase boundary for the latter is pre- 
dicted to occur at 1.49 T. The agreement between theory 
and experiment is very satisfying. 

Second, following the Clausius-Clapeyron equation 
(CC), dB m /dT = -AS/ AM, the transition at B m is 
connected with a strong entropy reduction of a thermody- 
namic phase transition subject to the value of dB m /dT. 
Because the transition at B m is hysteretic determina- 
tion of the phase boundary using CC depends on the 
choice of dB m /dT. As a first step it is instructive to 
suppose Kagome ice behaviour for B < B m . In this 
case a rigorous calculation predicts an entropy reduc- 
tion of 0.672 -1 mol _1 -ion [48]. The corresponding value 
of dB m /dT = 0.073 T/K, as shown by the solid black 
line passing through (7^,2?^) in Fig. [2] is perfectly con- 
sistent with our data and follows the cross-over line for 
T > T^. However, the phase boundary would lie asym- 
metrically with respect to 23+ and 23~. If we assume 
instead that B m is located midway in the hysteretic 
field range or at the field of the peak of dM/dB we 
find dB m /dT w 0.2 T/K. In this case the expected en- 
tropy release at B m exceeds the entropy associated with 
Kagome ice by an unphysically large factor of 3.65. Hence 
the latter prescription must be inappropriate. 

Strong support for the assumption of Kagome ice 
behaviour for B < B rn may be seen in the magne- 
tization, which forms a very well defined plateau at 
3.55 /^bHo^ 1 (cf Fig.|3](a)). This is quantitatively in 
excellent agreement with theoretical prediction. More- 
over, above B m the magnetization assumes a well-defined 
value of 5.20 /ib Ho -1 for T — > characteristic of a 
three-in/one-out spin configuration. Thus the magne- 
tization of Ho2Ti207 is quantitatively comparable with 
that of Dy2Ti2 07 and consistent with the dipolar spin 
ice model. The temperature dependence of 23+ and 2?~ 
and the magnetization provide unambiguous evidence of 
a first order transition with a large entropy reduction. 
Experimental determinations of the residual entropy of 
the Kagome ice state inferred from the specific heat in 
Dy 2 Ti 2 07 for B < B m are scattered between 0.44(8) - 
0.8(1) J K" 1 mol _1 -Dy PHI49], consistent with the value 
inferred from the magnetization, 0.5 JK _1 mol _1 -Dy [7]. 
The temperature dependence of the entropy release in 
the magnetocaloric effect in Dy2Ti2C>7 is attributed 
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FIG. 4: (Color online) Magnetic phase diagram of Ho2Ti207 
for magnetic field parallel to [111]. With respect to B = 
strong freezing and hysteresis emerges in the magnetization 
below Tg = 0.67 K in the presence of the residual dynamics 
seen microscopically [26H28] . At high fields a line of first 
order metamagnetic transitions terminates in a critical 
endpoint. The metamagnetic transition separates Kagome 
ice behaviour from a three-in/one-out configuration when 
going from below to above B m . 



to magnetic correlations [50] . 

We finally show in Fig. [4] the magnetic phase diagram 
of Ho2Ti 2 07 for the [111] axis inferred from our mag- 
netization data. In zero field Ho2Ti2C>7 enters a spin 
ice state with strong spin blocking below T *. A moder- 
ate field stabilises a magnetization plateau characteristic 
of Kagome ice. Further increasing the field results in a 
line of first order metamagnetic phase transitions up to 
a critical endpoint at 7 1 *, = 0.37 K and B'* n = 1.52 T. 
This line of transitions separates Kagome ice from the 
spin polarised 3-in/l-out configuration. While the criti- 
cal field for T — > scales experimentally with the J e ff, in 
agreement with theory, Tq and T* n appear to be material 
independent for Ho2Ti207 and Dy 2 Ti207. The tuning 
parameter for these energy scales remains to be explored. 

In conclusion, the remarkable analogy we observe be- 
tween the properties of Dy2Ti207 and Ho2Ti2C>7 es- 
tablishes the field induced liquid-gas like transition, 
interpreted in terms of the condensation of magnetic 
monopoles, as a more generic phenomenon within spin 
ice systems. However, given the importance of dipolar 
interactions and the resultant power law correlations |25j 
as an essential prerequisite for a description of the ex- 
citations in terms of magnetic monopoles, it also seems 
clear that further Ising like compounds not based on Ho 
or Dy, which have similar classical magnetic moments, 
must be investigated. In this way, the strength of the 
dipolar coupling could be controllably tuned. Perhaps 
most remarkably, the phase boundaries appear to be in- 
dependent of the strength of the hyperfine interactions, 
which are much stronger in Ho 2 Ti 2 07. It is important to 
note that the dc magnetization is a measure of the zero 
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frequency response of the system. Subleading transverse 
interactions J± which are material dependent and could 
lead for instance to differing quantum tunnelling ampli- 
tudes between spin-ice states are predicted to affect the 
finite frequency response [5T]. A complete description 
of the excitations in terms of magnetic monopoles must 
account for such differences. 
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SUPPLEMENTARY MATERIAL FOR: FIRST 
ORDER METAMAGNETIC TRANSITION IN 
HO2TI2O7 OBSERVED BY VIBRATING COIL 
MAGNETOMETRY AT MILLI-KELVIN 
TEMPERATURES 

In this supplement we present a detailed account of the 
precise temperature and field history used in our mea- 
surements and the terminology used to refer to specific 
temperature and field sweeps in the main text. 

Temperature dependence 

We distinguish two types of temperature sweeps, 
namely zero- field cooled/field heated (zfc/fh) and field- 
coolcd/field heated (fc/fh). In both cases data were 
recorded while heating the sample continuously at a rate 
of 5mK/min. Data denoted as zfc/fh were recorded af- 
ter initially cooling the superconducting magnet system 
and cryostat from room temperature, to guarantee the 
absence of any remanent magnetic fields. To record data 
denoted as fc/fh, we recooled the sample to base temper- 
ature after the zfc/fh measurement was completed (at 
roughly 1.5 K), keeping the magnetic field unchanged. 
Data were then recorded again while continuously heat- 
ing. 

Field dependence: Protocol A 

For field sweeps following protocol (A) the empty VCM 
was at first demagnetised and the sample cooled in zero 
magnetic field to T = 1.8 K. Data were subsequently 
recorded in three field sweeps as follows: (Al) — > +5 T, 
(A2) +5 T ->• -5 T, and (A3) -5 T -> +5 T. At the final 
field of +5 T the temperature was changed. At each sub- 
sequent temperature two field sweeps were carried out, 
namely (A4) +5T ->• -5T and (A5) -5T -> +5T. 



At temperatures above O.IK all data were recorded 
while sweeping the field continuously at 15mT/min. Be- 
low 0.1 K, to avoid heating effects due to eddy currents, 
data were recorded in a step mode, holding the magnetic 
field constant while recording the magnetisation. 

Field dependence: Protocol B 

For measurements following protocol (B) the empty 
VCM was at first demagnetised and the sample subse- 
quently cooled in zero magnetic field from room temper- 
ature directly to the temperature of interest. Data were 
then recorded in a sequence of three field sweeps: (Bl) 
-> +5T, (B2) +5T -> -5T, and (B3) -5T -> +5T. 
For the next set of measurements following protocol (B) 
the sample was removed from the VCM and the empty 
VCM again demagnetised. 

At temperature below 0.1 K and for all sweeps of type 
(Bl) data were recorded in a step mode. All other data 
were recorded while sweeping the field continuously at 
15mT/min. No differences were observed between (A4), 
(A5) and (B2), (B3), respectively. 





sweep # 


field values 


condition 


protocol (A) 


Al 


OT -»■ 5T 


zfc to 1.8 K 




A2 


+5T -> -5T 






A3 


-5T -> +5T 






A4 


+5T -> -5T 


field-cooled 




A5 


-5T -> +5T 


field-cooled 


protocol (B) 


Bl 


OT -> +5T 


zfc from ambient 




B2 


+5T -> -5T 






B3 


-5T -> +5T 





TABLE I: Summary of protocol (A) and (B) used in mea- 
surements as a function of magnetic field. Zero-field-cooled 
and field-cooled starting conditions are denoted as zfc and fc, 
respectively. 



